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Abstract 

In this paper we propose and analyze a feasible scheme where the detection of a single scattered 
photon from two trapped atoms or ions performs a conditional unitary operation on two qubits. 
As examples we consider the preparation of all four Bell states, the reverse operation that is a Bell 
measurement, and a CNOT gate. We study the effect of atomic motion and multiple scattering, 
by evaluating Bell inequalities violations, and by calculating the CNOT gate fidelity. 
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I. INTRODUCTION 



Implementing a quantum controlled-not (CNOT) gate is a key step in present attempts 
towards quantum computation [|l], 0]. Many different schemes for CNOT gates have been 
proposed J3], f|, ||, || [?[, and most of them require a strong quantum interaction between 
the particles that are used to carry the logical qubits. In practice, the quantum interaction 
is very often perturbed by classical noises, such as the individual motion of neutral atoms 
in a laser-induced potential well ||, or the collective motion of ions in a Paul trap |J|. 
Though order-of-magnitude estimations show that most of the CNOT gate schemes may be 
realized in principle, detailed analysis discovers many difficulties in eliminating all sources 
of classical noise for given experimental conditions. For instance, the perturbations due 
to thermal photons, photo-ionization, spontaneous emission..., make that the conditions 
for a fast CNOT gate operation through transient excitation to Rydberg states are only 
marginally satisfied || . If one looks at cavity- induced atom-atom coupling ( "cavity-assisted 
collisions" DTOf) in the optical domain, our estimations show that most schemes for cavity- 
enhanced coupling between the particles reliably works when g 2 /(hz , -f) > 10 3 , where g is 
a cavity mode- atom coupling constant, and k and 7 are respectively the cavity and the 
spontaneous emission damping rates. Though it is possible in principle to reach high value 
of g 2 / (/sry) fllTf , putting together very small high finesse cavities and reliable traps is far from 
straightforward. These considerations encourages us to look for "non-traditional" CNOT 
gate schemes, which do not require a direct interaction between particles, but rather use 
an interference effect and a measurement-induced state projection to create the desired 



operation [12fl . It was proposed in |L3[ |Tj] to create an entangled state of two atoms simply 
from the detection of a photon, spontaneously emitted by one of the atoms in such a way 
that the emitting atom can't be recognized. In such a scheme there is no direct interaction 
between atoms, and in principle the atoms can even be located very far from each other. 

In this paper we propose to extend the ideas of fl3f , fL4| , to realize a full quantum 
CNOT gate, or a Bell-state measurement, or more generally to implement conditional uni- 
tary operations. Our scheme will be based on an experimental setup using two atoms in 
two neighboring microscopic dipole traps || [15]], but it can be readily applied to other 



systems. In Section 1 we will describe how to realize a conditional unitary transformation 
that maps the four factorized states of two qubits onto the four maximally entangled Bell's 
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states. Since a convenient experimental signature of entanglement is the violation of Bell's 
inequalities (BI) |TB| , we will evaluate the result of a test of BI on the "transformed" pair 
of qubits, taking into account imperfections due to the motion of atoms (Section 2) and 
to the spontaneous emission of two photons by two atoms (Section 3). BI measurements 
are studied quantitatively in Section 4. In Section 5 we describe a CNOT gate based on 
the Bell's states created by the procedure of Section 1, and we calculate the fidelity of this 
gate, taking into account the motion of the atoms in the traps and the possible spontaneous 
emission of two photons. Finally we discuss these results and suggest developments of the 
proposed scheme. 

II. PREPARING FOUR ORTHOGONAL BELL'S STATES 

We consider two atoms i = 1, 2, trapped in two separate dipole traps, and prepared in 
one of two states |e) i or \g) i of the ground state hyperfine structure. We represent four 
initial states of the two-atom system as a vector-column 

{\a(3}} = {\gg) , \ge) , \eg) , |ee)}, a,(3 = e,g. (1) 

Each atom can be excited to one of upper states \e') i or \g') i by resonant a polarized laser 
fields of Rabi frequencies Q g i, VL eil as shown in Fig.l. The fields are weak, so that the 
probability to excite both atoms is much smaller than the probability to excite only one 
atom. An excited atom may emit spontaneously a photon, with the wave vector k and 
certain polarization, on 7r-polarized \g') i — > \g) i or |e'). — > \e) i transitions. Occasionally a 
photon passes through the optical system shown in Fig. 2 and it is registered by the photo- 
detector. We assume that the polarizer P transmits only 7r-polarized photons, and thus 
a polarized photons emitted on the \e') i — > \g) i and \g') i — > \e) i transitions will not be 
registered. 

After the excitation, the wave function of two atoms is changed from \ot(3) to \^ a )i ^13)2 
l*e>, = \e) t + b \g') i e 'CW*-HP->, W< = |^ + b | e >). e (Mr i+< M ^. = j^. + ^ (2) 

where 5ri describes fluctuations in the position of atom i near the equilibrium due to the 
motion of the atom in the trap, r, is the atom position at equilibrium, k e 5 are the wave 
vectors of the laser field resonant to either e — > g' or g — > e' transitions, tp° ai is the phase of 
the laser field fl ai , b 1 is a real constant. 
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FIG. 1: Scheme of the relevant atom transitions 




FIG. 2: Proposed scheme for conditional quantum logic. Atoms 1, 2 are placed in a focal plane 
of the input lens L of the optical system. The atom pair is excited by the laser field with the wave 
vector ki, circularly polarized in the yz-plane, and emits a photon with the wave vector k on the 
x-polarized transition. The polarizer P selects x-polarized photons, that are transmitted through 
an interferometer / towards a photo-detector. A mirror m of the interferometer is tilted, so that 
the images 1', 2' of the two atoms coincide on the photo-detector. The unit vector in the direction 
of a 7r-polarized atomic dipole is denoted as d, and 9q is the aperture angle of the lens L. 

The registration of a photon means that the wave function |\P a ) |\l>g) is projected to a Bell 
state \B a p) = |k) B \af3) where |k) is the state of the field with one spontaneously emitted 
photon. For example, \^ g ) \*5> g ) is projected to 

|k> B \gg) = |k) \B gg ) , \B gg ) = (1/V2) {\ge} e %^+^( k )+^l + \ eg ) e i[q,*r 1 +fci 1 (k)+ Vfll ] j ? 

(3) 

where q Q = k a — k, k = |k| and /j(k) is the optical length which a photon travels through 
the optical system towards the photo-detector. The optical system is set in such a way, 
that images 1', 2' of atoms 1 and 2 perfectly coincide on the photo-detector. This means 
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that Zj(k) is the same for all registered photons, and therefore feZj(k) = kk. Introducing the 
vector-column {\B a p)} = {\B gg ) , \B ge ) , \B eg ) , \B ee )} of the Bell states we can express them 
in terms of the initial states ([!]) as {\B a p)} = [B]'{\a(3)}, where 



[B]' 



V2 



Q e «(q 9 '5r2+</3 g 2+fc'2) e i(qgSr 1 +ipg 1 +kl 1 ) g 

e i(q e (5r2+(/Pe2+fc^2) Q g gi(q 9 <5ri+</3 9 i+fcii) 

gi(q e 5ri+v?el+Wl) g g gi(q 9 <5r2+y 9 2+fc£2) 

g gi(q e <5ri+</j e i+fcZi) g«(q e <5r 2 +</3 e 2+fci2) g 



(4) 



is a matrix of Bell operator B. 

In general, the wave function \B gg ) (\B ge )) is not orthogonal to \B ee ) (\B eg )). In order to 
make sure that all Bell's states are orthogonal one has to satisfy two conditions 







(B ee \B gg ) = W ge e *(v»i-w+«i-«2) + W * g e 



i(<Pg2-<Pel+kl2 — kh) 



= (B eg \B ge ) = W gg e^-^+kh-kh) + ^ gi(c, e2 - Ve i+ 



kl2—kli) 



(5) 
(6) 



where W ai g = e *( qa<5ri ~ q ' 3 ' 5r2 ). If the atoms are very cold in a steep trap, so that they are 
deeply in the Lamb-Dicke regime, one should take W a p = 1 in Eqs.@ and (^). But we point 
out that the resulting conditions are actually independent of the atoms motion. Indeed, 



<Wap) 



w: 



V a V p = l-D(T), V a 



\n=o (2n) 



2n\ 



(7) 



q a 



where (...) T and (•••) qQ means, respectively, the average over the atom motion and over 
directions of registered photons. We average separately over symmetrical and statistically 
independent motion of each atom, drop index % in Sr^, and introduce parameter D(T), 
< D(T) < 1, where T is the temperature associated with the random motion of the 
atoms. One can see that {W a p) and (w*^ disappear from orthogonality conditions (|J), 
(0), which are reduced to a single condition 



<P% ~ <Pgl + &2 ~ Vll + (ke + k 9 )(r 2 - n) + 2k{l 2 - h) = 7T 

In our geometry we have (p° al = ip° a2) and thus this condition becomes 



(k e + k 9 )(r 2 -r 1 ) + 2A;(/ 2 -/ 1 ) = tt. 



(8) 



(9) 



There are various ways to fulfill this condition. If one chooses k e = — k 9 , i.e., the cr+ and a_ 
lasers are propagating in opposite directions, the condition for orthogonality of the Bell's 



states is obtained by adjusting the interferometer path difference so that k(l 2 — h) = tt/2. 
But it is also possible to take k e = k g = k^, together with k^(r 2 — ri) = ir/2, obtained by 
adjusting the trap's positions. Assuming then that kl 2 = kl\ = 2mr for simplicity, taking 
(here and everywhere below) the origin of the coordinate system on the atom 1 and defining 
q = k^ — k, the Bell operator matrix can be written: 

ie^ 5r2 e iq5ri 



W = T2 



iq<5r2 



3 iq<5ri 













3 iq5ri 



ie 



jq<5r2 



e lcldri ie lcldr2 

which converts four initial atom states to four Bell states, which are orthogonal in average 
over the atom motion. Though the condition for the orthogonality in average depends only on 
the atoms equilibrium positions, the final fidelity of the conditional unitary transformation 
will obviously depend on the atoms motion, due to the Sri and 5r 2 in the [B]' matrix. 

In order to simplify the local operations used in the rest of the paper, it is convenient 
to perform two phase transformations for atom 2, that make the change |e) 2 — > — % \e) 2 just 
before the photon observation and |e) 2 — > i\e) 2 right after it. Taking into account such 
transfornations as diag{l, — i}[B]'diag{l, i, 1, i} = [B], where diag means diagonal 
matrix, we find 



[B] 



V2 



Q _ e *q<5r2 e «q<5ri g 



3 iq<5r 2 



3 iq<5ri 










3 jq<5ri 






iq5r2 



3 iq<5ri 



_ e iq<5r 2 







(10) 



which has real elements in the absence of atom motion 5ri = 0. In a geometry where the 
phases (p° ai can be independantly controlled, one can obtain the matrix (III) more straight- 
forwardly, for example by choosing in Eq.([8D kl 2 = kl\ = 2mr, k e = k 9 = k^, ri = 0, 
and 

¥^ = ¥^=0 ^ 2 = -k L r 2 , ^ 2 = 7r-k L r 2 . (11) 

Below we refer to [B] as a Bell operator matrix supposing either that condition (|]) is true and 
the Bell operation is the photon observation procedure with the two phase transformations 
for atom 2, or that there is only the photon observation, but conditions (JTTJ) are satisfied. 

In order to get a physical understanding about the quality of the Bell states preparation, 
we will now look in detail whether the prepared states can violate Bell's inequalities. In 



these calculations we will use the expression 
similar results could be easily obtained by in 
matched situation where the atoms' positions 
experimental geometry). 



fllOD corresponding to k e = k g = k^, but 
the case where k e = — k 5 (the fully phase- 
would cancel out is not accessible with our 



III. BELL'S INEQUALITIES 

When the atoms are prepared in a Bell state, the statistical behavior of measurable 
quantities (such as the population of state is governed by the entangled wave-function 
\B a p). Here BI will be used as a simple experimental characterization of the degree of 
entanglement of the atom pair. As we will see below, either atoms motion or simultaneous 
excitation of the two atoms may reduce or even suppress the BI violation. 

In order to test BI we carry out the following sequence of operations 

1. The atoms are prepared in one of states ([!]); 

2. Atoms are excited by a weak laser pulses under the conditions of Eq.(^); 

3. One spontaneously emitted photon is registered. If there is no photon after some 
delay, the stages 1, 2 are repeated until one photon is registered; 

4. Raman transitions for each atom are carried out so that 

\g) i -+caa(6 i )\g) i -w\(6 i )\e) i \e) i ->• cos (0 4 ) |e) i + sin (0<) \g) i ; (12) 

5. Populations of \a) i states are measured; 

6. Operations 1-5 have to be repeated until a full statistical ensemble of results for the 
population of \a) i states is obtained; 

7. The steps 1-6 are repeated for four different Raman transitions with four pairs of 
angles {9 1 ,9 2 }, {9[,9' 2 }, {9[,9 2 }, and {9 1 ,9' 2 }. 

After the operations 1-6 are carried out, the state of atoms and a photon is |k) RB \af3), 
where the operator R describes Raman transitions ( W% for two atoms. The matrix [RB] 
of RB operator is the matrix product [-B][i?(#j)] where the matrix \R(9i)] for the Raman 
transitions is given by Eq.(44) of Appendix 2, and [B] is given by Eq . (|l0|) . Here and below 



we denote the dependence on 9\ and # 2 as a dependence on 9i, when it does not lead to 
confusions. 



7 



Let us call P 



(75), 



a g \yi) the probability to find atoms in state [jS), while the initial atom 
state is \a/3). By taking the modulus square of each matrix element in [RB] one can find 



(ee). 



99 



0.5 



sin 



2 ) + 0.5Q(T, 0i 



P^(9 i ) = Pjf >(*«) = P^(9 i ) = P^m = 0.5 [cos 2 (0 a - 9 2 ) - 0.5Q(T, 9 t 

Pifm = pt ] m = Pifm = p£°(0<) = 0.5 U 2 (9 1 + e 2 ) + o.sq(t, 9, 



(13) 



(ee), 
ee < 



0.5 



sin 2 (0x + 2 ) - 0.5Q(T, 9 { 



eg \"U - eg 

where Q(T, 0*) = £>(T) sin (20 x ) sin (20 2 ) and D(T) is given by Eq.(0). One has D(T) = 
in the absence of atoms motion, in this case the maximum violation of Bell's inequalities is 
obtained. In the opposite (high temperature) situation, where D{T) — > 1, one obtains from 
Eqs.( 



P { gl 9) m 



0.5 



sin 2 (0x) cos 2 (0 2 ) 



sin 2 (0 2 ) cos 2 (0i ) 



(14) 



and similar expressions for the other probabilities. This corresponds to a "classical" limit 
where BI cannot be violated. Thus D(T) is a "decoherence parameter" which grows up with 
the temperature from D(0) =0 to maxD = 1. 

For each atom i = 1, 2 we define a random variable £j, with values +1 or —1 depending 
on whether an atom is found, respectively, in \g) i or \e) i state after registering a photon and 
carrying out the Raman transition. With the help of Eqs.fll^) one can find = (£2) = 0, 
(£1) = (£2) = 1) where the average is made over the results of a sequence of operations 1 - 
6. The correlation functions are given by: 

(£i£ 2 ) - (£1) (£2) 



E a j3{9i) = 



(£i£ 2 ) = 2[P* 



% e) ( 



so that 



E ge (0i) = -E gg {6i) 
E eg (9i) = —E ee (9i) 
As usual we define the quantity 



SafiiQi, = E, 



cos [2(9 1 -9 2 ))-Q(T,9 i ) 
= cos [2(9 1 + 9 2 )]+Q(T,9 i ). 

E a p(0i9' 2 ) + E a/3 (9' 1 9 2 ) + E a j 3 (9 i ) 



(15) 



(16) 



for each initial state \a,/3), then the Bell's inequalities read [17] 



-2<s af) (e i ,e' i )<2. 



(17) 



The results (|T3|) - (|T6|) are similar to ones obtained for the BI test with polarization-entangled 
photon pairs [|T^] , the difference is that here the decoherence is taken into account by means 
of D(T). In the next Section we look for the violation of inequalities ( |TTD for each initial 
state of the two-atom system. 

IV. EFFECT OF ATOM MOTION ON BELL'S INEQUALITIES TEST 

In order to predict accurately the value of SapiQi, 9^) we have to calculate the factor D(T), 
which depends on the trapping potential and the aperture angle 9 of the input lens of the 
optical system. In general, the trapping potential is an-harmonic, non-symmetric and 9 is 
not small. All of these complicates the precise calculation of D(T), which will be carried out 
elsewhere. Here we will consider a simple order-of-magnitude estimation, using an harmonic 
approximation for the trapping potential. The procedure carried out in Appendix 1 (see 
also |Lj|) leads to 

D{T)Kl-e- T ' T -, k B T cr = ^l, (18) 

where T cr is a critical temperature such that D(T > T CT ) « 1, v e ff is an effective frequency 
of atom motion in the trap. For the aperture angle 9q = tt/4, as it is in our case, v~h ~ 
1.25^J 2 + 0.75z/||~ 2 , where u± and v\\ are the frequencies of the motion of atoms in x, y 
and in z directions, respectively; Er = h 2 k 2 /(2m a t) = hvR is the recoil energy, and ks is 
Boltzmann constant. Using vr = 3.6kHz for Rb 87 atoms and our estimations u± = 200 kHz 
and = 50 kHz we obtain u e ff = 55 kHz and T cr « 20 /iK. In general, u e ff and T cr depend 
on #o and the direction of the laser field. The maximum T cr is reached when k is parallel 
to ki for the most of the emitted photons. For the geometrical arrangement displayed in 
Fig. 2, the variation of T cr as a function of 9 Q is given by Eqs. (|42]) of Appendix 1, and it is 
displayed in Fig. 3. 

Let us choose parameters of Raman transitions 

Q x = 0, 9 2 = x, 9\ = 2x, 9' 2 = 3x, (19) 

and T/T cr = 0.5, for such case factors S ge (x) = —S ee (x) and S eg (x) = —S gg (x) are shown 
in Fig. 4a. One can observe the violation of BI \S a/ 3(x)\ > 2 for S ge (x) and S ee (x), while 
BI are satisfied for S eg (x) and S gg (x). This situation can be inverted by choosing 6>i = 0, 
6> 2 = — x, 9[ = 2x, and 9' 2 = —3x, so that BI will be violated for S eg (x) and S gg (x) but 
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FIG. 3: Critical temperature as a function of the aperture angle of the optical system. 




FIG. 4: (a): S ge ( x) — — S ee {x") (solid line) and S eg (x) — — S gg (x) (dashed line) at T /T cr — 0.5 
and angles Q\, 62, 0[, 0' 2 given by Eq.(|l^). Bell's inequalities are violated for S ge (x) = —S ee (x). 
(b) The maxima of \S a p\ for different T/T cr . Dotted lines in Fig.4a mark \S\ = 2 and |»S| = 2\/2. 

satisfied for S ge (x) and S ee (x). Therefore all four states do violate BI, but the combination 
of angles to be used depend on the state in the pairwise fashion just described. Fig. 4b shows 
the maxima of \S a p\ versus the normalized temperature of the atom motion found for # 12 , 
8[ 2 given by Eqs.(|l9D. The condition to violate BI for all four states (for suitable choices of 
Raman angles) is therefore that T < T cr . 

V. EFFECT OF MULTIPLE SCATTERING ON BELL'S INEQUALITIES TEST. 

Now we take into account the excitation of two atoms together and examine how it 
influences the BI violation. Let us first consider the case of only three levels in each atom 
shown in Fig.5. We examine a possibility that a photon emitted by one atom is registered, 
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FIG. 5: The event when two atoms are excited together, both of them came to the |e) state, but 
one emitted photon is missed. 

while another atom also emits a photon on |e') — > |e) transition, but this photon is missed. 
After the excitation to \e') i state the wave function of atom i — 1, 2 is given by Eqs.(||) 



\^g) 1 = a\g) 1 + b\ef) 1 e 



™ r \ \y g ) 2 = a\g) 2 -b\e') 2 e^, 



where we suppose that conditions ([TT|) are satisfied. If the photon k is emitted by one atom 
and registered, while another photon k' is emitted by the other atom and missed, then the 
atoms go from \e'e') to |ee) state and we have: 

f = e *(q<5ri+q'<5r 2 +v?2) _|_ e «(q<5r2+q'5ri+</j' 1 ) ^0) 

where q' = k^ — k', ip[ is the phase of a missed photon emitted by atom i — 1,2 and we 
suppose, as usual, fc/i = kl 2 = 2irn. It may also happen that one atom emits the missed 
photon first, and then the registered photon comes from another atom. In that case, one 
has to change the field state |k, k') in Eq. (f20|) to |k',k). After registering a photon and 
performing the phase transformation |e) 2 — > i |e) 2 , the state \gg) is projected to the state 

1 r 'jrtn {eg) _ ^5r 2 lge) \ |k) + (0 l/ 2/ | ee) (|k / ; k) + |k; k , }) 



(21) 

where £ = (b/a) 2 . Taking into account that the field state |k', k) is orthogonal to |k, k'), 
and calculating (|/| 2 ) T = 2, one obtains the normalizing factor iV = 2[1 + 2£]. 

After carrying out the Raman transitions, the atom states in the right part of Eq.(^Il) 
are changed in accordance with the transformation (|T2"|) . Following the procedure of Section 
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Ill we find 



E, 



-[1 -£>(T)]sin20i sin2# 2 



cos 28i cos 26*2 
1 + 2£ 



(22) 



Fig. 6 shows S gg (x) calculated with the help of Eqs.(|l6l),(p2|) for T/T cr = 0.5, 9i, 2 , 9[ and 
9' 2 given by Eqs. (|l9D and various £. If the state \e') i is excited by a weak "square" pulse, so 
that Q g i is constant during the excitation time and zero otherwise, then £ = |O g j| 2 /<5 2 <C 1, 
where 5 is the detuning from the resonance on \g) i — > \e') i transition. According with Fig.6, 
BI are still violated for \tt gi \ 2 /5 2 < 0.15. 




0.2 0.4 0.6 



x/k 



FIG. 6: Factor S gg {x) for Raman transition parameters given by Eq.(19), T/T cr = 0.5 and £ = 
(curve 1), £ = 0.05 (2), £ = 0.15 (3) and £ = 1 (4). BI violation is possible for £ < 0.15. 



It is convenient to write the final two-atom state, taking into account simultaneous exci- 
tation of two atoms for each initial state \af3), under the following matrix form: 



2/N[\k)B + |k,k')S (2) ) \a{3) 



(23) 



where the matrix [B] of the operator B is given by Eq. (PH|) , the matrix of the operator B^ 
is 

r o 1 
10 
10 
10 

and the state |k, k') is orthogonal to |k) and normalized to 1. 



[B 



(2)1 



(24) 
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VI. BELL'S STATE MEASUREMENT 



In the previous section we have shown that four orthogonal Bell's states can be prepared 
from four initial factorized states, under the condition of detecting a single photon. The 
reverse process, usually known as a Bell measurement, is actually also realized using the same 
scheme. Here we will evaluate the efficiency of a whole sequence, including the preparation 
followed by the measurement - two successive clicks will be therefore required. 

By carrying out the steps 1 - 3 of the procedure described in Section |T| we prepare a 
Bell state \B a/3 ) of two atoms. Then \B a/3 ) can be projected to the pure state \aj3) of two 
atoms, or "measured" , by proceeding the steps 1-3 with the phases of the laser fields 

<p° e2 = 7T - k L Ar, <p° gl = V ° el = 0, <^ 2 = -k L Ar. 

Taking into account the multiple scattering, one arrives to the final state of two atoms and 
spontaneously emitted photons after the Bell's state preparation from \af3) state followed 
by the Bell's state measurement 



sjl/fi ( k) h + k, k') B {2) ) (|k> B + |k, k') 5 (2) ) \a(3) , (25) 



where N is normalizing factor. A matrix [B] of an operator B is [B]^ 1 with Sri replaced by 
8i i, operator B^ is the same for the Bell's state preparation and the measurement. Because 
of the preparation and the measurement of a Bell's state are separated in time, all field states 
in Eq . (|25|) are orthogonal to each other and the average over the atom motion (<5rj<5fj) = 0, 
i,j = 1,2, since the atom motions on different time intervals are not correlated. 

The Bell's state measurement is not perfect due to the atom motion and the multiple 
scattering, so that the state fl25|) is, in general, a linear combination of four states (|1|). If a 
fidelity of the Bell's state measurement is high, the probability to find atoms in \a(3) initial 
state after the measurement approaches 1, while the probabilities to find any other atom 
states tends to 0. A matrix for the transformation of the vector-column of states ([I]) after 
the Bell's state preparation and the measurement is 



1/N {[B] [B] + [B®] [B] + [B] [B^] + [B^f ) . (26) 

Taking the square modulus of each element in the matrix ( p6|) and calculating N = [1 + 2£] 2 
one converts matrix ( p6|) to a matrix of probabilities to find atoms in \^yS) final state starting 
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with \a(3) initial state 



;i + 2iy 



(27) 



l-h+Ae 2£ 2£ f x 

2£ A i-A+4e 2 2£ 

fi 2£ 2£ i-fi+^e 

f 1 (T) = D(T)-D 2 (T)/2. 

In the case of perfect Bell's state preparation and the measurement the matrix fl2"T|) has 
diagonal elements equal to 1 and other elements equal to 0. Thus, we can take the diagonal 
element of the matrix (|27|) as the fidelity Fb(£,T) of Bell's state measurement 

A£ 2 + D(T) -D 2 (T)/2 



F B (£,T) = 1- 



(28) 



(1 + 20 2 

Fidelity Fb is shown in Fig. 7a as a function of T/T cr for various £, it is shown in Fig. 7b 
as a function of £ for various T/T cr . Small increase in Fb for large £ « 1 is because of the 




T/T,.,. 




FIG. 7: (a) Fidelity factor Fb for Bell's state measurement as a function of T/T cr for £ = (curve 
1), f = 0.05 (2), £ = 0.15 (3) and f = 1 (4); (b) F as a function of £ at T/T cr = (curve 1), 
T/T cr = 0.2 (2), T/T cr = 0.5 (3), T/T cr = 1 (4). 

contribution of processes \aa) — > \/3/3) — > |aa), a ^ (3 grows up with £ due to the multiple 
photon scattering. However this is not so important for practical cases, where £ <C 1. 



VII. QUANTUM CNOT GATE. 

We have shown so far that conditional Bell states preparation and measurement can 
be successfully achieved. Our result is actually more general than that, and shows that 
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arbitrary "conditional" unitary transformations on two qubits can be achieved by using 
Raman rotations (applied locally to each atom) and the detection of a single click. In order 
to demonstrate this we will now show that the four Bell's states preparation \B a p) = B \af3) 
can be turned into a "controlled- not" (CNOT) operation C, described by the matrix 



[C] 



10 
10 
1 
10 



(29) 



We prove that in our case C = H 2 BHi, where Hi 2 are some local (single-atom) operations 
and the matrix of Bell operation B is given by Eq. (|10"D . The CNOT operation can thus be 
realized by the following procedure: 

1. One of initial states ([I]) of atoms is prepared; 

2. The local operation Hi is carried out; 

3. Atoms are excited and a spontaneously emitted photon is registered. If there is no 
photons registered for a time t 3> operations 1-3 has to be repeated; 

4. The local operation H 2 is carried out. 

Let us suppose, for a while, that atoms do not move, so that the Bell operator is Bq, which 
matrix [Bo] is given by Eq.(TT0|) with 5ri = 0. By definition [C] = [Hi][B ][H2], where [#1,2] 
are the matrices of local operations Hi 2 , and therefore 



[H 1 ) = [C\[H 2 )- 1 [B 



(30) 



Taking the matrix [-f^] -1 as a general local transformation for two-level atom, inserting it 
in Eq.(^) with the requirement that the matrix product on the right of Eq. ([!(]) should be 
a local transformation, we obtain 



— 1 li — 

iii i 
-11-1 1 








-1 


-i 


1 


i 





1 


72 





-1 


-i 




—% 





1 



(31) 



Details of the procedure of determining of [#1,2] are given in Appendix 2. As it can be seen 
from Eqs. (|5T|) , (|52f) of Appendix 2, operation Hi is the phase transformation \g) 2 — > i \g) 2 , 



after which the Raman transition (|i~2|) with 61 = 7r/4, 8 2 = —n/A is carried out. Operation 
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H.2 starts with the Raman transition (|T^) with 9\ = —tt/4:, 8 2 = —ti/2 after which one makes 
the phase transformations |e} 1 — > —i \e) v \g) 2 — > —i \g) 2 . 

Now we take into account the atom motion, the simultaneous excitation of two atoms and 
find a fidelity of CNOT operation. We suppose, that transformations are much faster 
than a period of the atom motion in the trap, in such case B.\^ does not depend at all on 
the atom motion. Indeed, by carrying out a fast local operation with atom i, one can chose 
the origin of the coordinate system in that atom, which means r i0 + 5r = 0. Thus, using 
formula flltlD with Sr^ ^ and Eq.(|23|) we obtain an operator C of a non-perfect CNOT 
transformation 

' (32) 



C = ^2/N (|k> C {1) + |k, k') C {2) ) , 
where matrices of operators C**- 1 ' 2 ^ are 



[CW] = m ■ [B] ■ [B 



2J- 



[C^] = [H,] ■ [B^] ■ [H 2 \. 



(33) 



1, 



The matrices are given by Eqs.flSID, and matrices [B], [B 2 ] are given by Eqs.(|10|) 
respectively. 

We can build now matrices [|C^' 2 ^| 2 ], which elements are the square modulus of respective 
elements of [C^ 1 ' 2 -'], and calculate the matrix [C p ] of probabilities to find atoms in \^y5) state 
after CNOT operation, while \a(3) was the initial atom state 



2 

N 



{[|C«| 2 ] + [|C«| 2 ]} 



1 + F 1-F 

1 - F 1 + F 

1-F1+F 

1 + F 1-F 



where 



D(T) 



1 + 2^ 



(34) 



< F < 1 is the fidelity of CNOT operation d32|). Factor F is shown in Fig.8a as a function 
of T /T cr for various £, it is shown in Fig.8b as a function of £ for various T/T cr . 



VIII. DISCUSSION 



An important problem in the experimental demonstration of the conditional quantum 
CNOT gate operation is the suppression of the atom motion, which can be done by cooling 
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FIG. 8: (a) Fidelity factor F for quantum CNOT gate as a function of T/T cr for £ = (curve 1), 
£ = 0.05 (2), £ = 0.1 (3) and £ = 1 (4); (b) F as a function of £ at T/T cr = (curve 1), T/T cr = 0.2 
(2), T/T cr = 0.5 (3), T/T cr = 1 (4). 

atoms in the traps up to temperatures of few /zK, or by increasing the atom oscillation 
frequencies. The last can be done, for example, by using standing-wave trapping fields, that 
separates a trapped potential into several narrow wells with oscillation frequencies much 
higher than the present 50 kHz obtained for the longitudinal motion in a tightly focussed 
beam. Multiple scattering gives rather small contribution, even for few percents of the atom 
excitation. An important characteristics of the fidelity of the CNOT gate operation is the 
decoherence parameter D(T). For a reliable theoretical determination of D(T), one needs to 
know accurately the trapping potential. However, before doing experiments on BI violation, 
D(T) can also be determined experimentally by looking at the interference fringes on the 
light emitted by the two atoms, irradiated on a closed transition [19| . 



Neutral atoms in a dipole trap are not the only candidates for implementing our condi- 
tional CNOT gate. In principle, the gate can be realized with other resonant objects such as, 
for example, trapped ions, or quantum dot molecules (QDM) incorporated in a solid matrix 
p0| , [21]. Each QDM consists in two closely positioned quantum dots with the ground state 
of each dot split into two or more close states. The advantages of QDM are their fixed 
positions in the matrix, and the possiblity to prepare the initial states electronically. The 
difficulty, however, is in providing the coherence during the gate operation, which is quickly 
destroyed by the electron-phonon interaction. 

The proposed simple scheme can be generalized straightforwardly to more complicated 
schemes with many elementary gates, which may be called "integrated conditional quantum 
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logic blocks", or ICQLB. They can be constructed by the increase of the number of atoms 
and (or) the number of ground states available in a single atom. There are p n initial states 
of n atoms, if an identical photon can be emitted on transitions to p different ground states. 
However, because of the photon observation process is not Hermitian, the maximum number 
N(n,p) of obtained orthogonal Bell's states is, in general, less than p n , though N(n,p) 
increases with n and p. One can see, for example, that only 7 orthogonal Bell's states are 
possible for three atoms with the level scheme of Fig.l, for any choice of phases of laser 
fields. Determination of N(n,p) is an important question for the theoretical modeling of 
ICQLB. 

Another theoretical problem is how to find local transformations, which convert N(n,p) x 
N(n,p) matrix of the generalized Bell transformation, obtained by the photon observation 
procedure, to the matrix of desirable logic transformation. The procedure of Appendix 2 can 
be generalized, in principle, to higher-dimensional cases, however it seems too cumbersome, 
and the development of a simpler procedure would be quite helpful. We underline that even 
complicated ICQLB operates in the same five steps as the CNOT gate described above. The 
step 1 is the preparation of initial N(n,p) states of atoms; 2 includes local transformations; 
3 is the excitation of atoms by weak resonant fields, repeated may be several times until a 
spontaneously emitted photon is registered; 4 is another the local transformation; 5 is the 
determination of final populations of atomic states. Each step can be carried out simultane- 
ously for all atoms together, so that the operation time of ICQLB is not much longer than 
for the elementary CNOT gate. The increase in the operation time of complicated ICQLB 
can happen however, because the probability for n atoms to emit more that one photon 
increases with n, so that lower intensities of the exciting fields are required in order to avoid 
multiple scattering. 

As a conclusion, though we do not propose here a specific way to make the present scheme 
scalable, it would be very interesting to study up to which point quantum computations may 
be realized using conditional logical elements such as the ones described above. 
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Appendix 1 



We consider an atom in the trapping potential as a harmonic oscillator, with a deviation 
from the equilibrium position given by 5vi = {5r x ,5r y ,5r z }. Then ^e* q<5r,; ^ = e~^ Sv ^ )t/ 2 ; 
which is the consequence that the thermal fluctuations of the position of a harmonic oscillator 
are described by the Gaussian distribution function, and therefore 



1 - D(T) = (e 



,iq(<5ri-<5r 2 ) 



-((q<5r) 2 ) 3 



(35) 



In the coordinate system shown in Fig. 2 k x = ksin9 cos (p^ ky — Jv sin#sin<y9, k z = kcos9 
and k Lx « k, k^ y = k Lz = 0, so that 

(36) 



[1 — sin 9 cosy?) (Sr x J +sin ^sin <p (5r y J + cos 9{5r z/ 



For the one- dimensional quantum harmonic oscillator with the mass m at , which oscillates 
along the axes £ with the frequency u^, £ = x, y, z in the thermal equilibrium 



8r\ 



coth 



Hue 



k B T 



't 2m at uj^ \2k B T ) m at ujf 
at Tioj^I '(k B T) <C 1. Taking u x = u y = 2ttv± and uj z = 2ttvu we obtain 

2 k B T 



(37) 



The average 



(q5r) 2 ) =t 

/ r 27T 2 m at 



(1 + sin 2 9 — 2 sin cos u?) —5- + cos 2 9^ 







2tt r o 

dip / sin9d9V(9,(p)u(9,(f), 
Jo 



where u(9, (p) is some function and 



V{9,(p) = C (l - sin 2 9 cos 2 ip) 



(38) 



(39) 



(40) 



is a probability that a photon, emitted to the solid angle of the optical system has a direction 
k, Co is the normalizing constant 



— = / dip sin 9d9{l — sin 2 9 cos 2 ip) 
Eq.(HDp is obtained from the relation 



47T 

T 



1 - -(3 cos O + cos 3 o ) 



(41) 



P(0 )¥ >) = Co £ [d ■ ^(k/A;)] 2 = C [l - (d ■ k) 2 /fc 2 ], 



A=l,2 
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where e,\(k) is the unit wave vector along one of two possible polarizations A = 1,2 of a 
photon and d is a unit wave vector along the direction of the dipole momentum of atom 
transition directed alone shown in Fie. 2. 

Thus, in order to find D(T) one has to insert Eq.(^) into Eq.fl3"5|) and calculate (...) q as 
shown by Eq.(]3"9"D. For a high input aperture of the optical system, that is our case, this 
procedure can hardly lead to an analytical result. For the estimations we use an approxi- 
mation 



■ 



which is as better, as D(T) is smaller. Thus, we arrive to D{T) = 1 — e T / Tar ( e °) where 

kT (0)- h2v "ff^ 1 A \\%) , MOo) U2] 



2E R ' vl ff {6 Q ) v\ v\ 



A± (6 ) = 1 I 47rC °( g °) U _ 5cosfl -cos 5 ^ A ^ = 87rC (6 ) ^ _ 5 cos 3 6 + 3 cos 5 6 * 
and C (9 ) is determined by Eq. fl41|) . 



Appendix 2 

General local transformations, which mixes the states and \e) i of the two-level atom 
i = 1,2, can be written in the form 

\g)i -> \g)i - Sie** \e) t )e 1 ^ |e>. -> (ae** \e) i + Sl e 1 ^ \g) t )e 1 ^. (43) 

They consist the phase transformation which changes the phase of the state \a) i — > e l ^ ai \ct) i: 
the Raman transition given by Eq . fll2"D with q = cos(6 l j), Sj = sin(^) and another phase 
transformation \a) i —>■ e l ^ ai For our purposes, however, it is enough to consider trans- 
formations (|43"D (or reverse to them) with <ft a i = 0. The matrix of the local transformations 
carried out with two atoms is, therefore 

[L{6U«d] = [M{U)Wi)\ = Mfa)) 

where 

[M(fa)] = [M{£ gl , U, e 3 2, U)} = diag{e*»,e*>% e*-»,e*"}, fa = fa + ^ 2 , (45) 



CiC 2 — CiS 2 — S\C 2 S\S 2 

ClS 2 CiC 2 — SiS 2 — S\C2 

s l c 2 ~ s l s 2 c l c 2 — c l s 2 

S\S2 S\C2 C\S2 CiC 2 



(44) 
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is the diagonal matrix of the phase transformation and [R(8i)} = [R{0i, 82)} is the matrix of 
the Raman transformation ([12]). 

Let us take [i^] -1 = [L(9i,^ ai )}, insert it into Eq.(pO[) and find 



[Hi] = [M^ai)} 



CiS 2 — SiC 2 CiC 2 + SiS 2 CiC 2 — S1S2 — CiS 2 — S1C2 

-CiC 2 — S1S2 CiS 2 — SiC 2 CiS 2 + SiC 2 CiC 2 — SiS 2 

CiS 2 — SiC 2 SiS 2 + CiC 2 SiS 2 — C1C2 S1C2 + CiS 2 

S1S2 + C1C2 S\C2 — C\S2 S\C2 + C\S2 C\C2 — S\S2 



(46) 



where [Mi(£ a j)] is a diagonal matrix obtained by interchanging two last elements of [M(£ Q j)]. 

Our goal is to determine 6i and £ a j, such that the matrix given by Eq.(f46"D can be repre- 
sented as 

[H 1 ] = [L(9 l ,U)] (47) 

with some 9 12 , £, a i- By comparing the matrices given by Eq. (pH|) and ( pE6] ) we see, that 
Eq.(f47"D can be true only if |cj| = \si\, that is when 8i = ±7r/4, 9i = ±37r/4; or when c.j = 



or Si = 0, while \cj\ = \sj\, j 7^ i. We chose C2 = 0, S2 = 1, that is for 6 2 = tt/2 and c\ = s± 
that is for 61 = 7r/4, so that 



[Hi 



p*£»e p^ge pi£ge p i£,g 



g^£ee g^£ee g^£ee g^£ee 
pi^eg pi^eg p^eg p^e 



(48) 



We can see now, that the matrix [L(#j,£ Q j)] is very similar to the matrix 
6*i = 7r/4 and #2 = — tt/4, so that 



HI) if we take 



[L(7r/4,-7r/4,e m )] 



3*?99 p*?99 p*?9 



3*?9 



g*?fle g*?ffe g*?9e g*?ge 

g^£eg p^eg pi£eg p^^eg 
g^£ee gi£ee gi£ee g^ ee 



(49) 



with 

4/3 = fal + 1/32- (50) 

Apart of the notations for phases, the only difference between the matrices in Eqs.([49"D 
and Eqs.(^HJ) is the opposite signs of the elements in the last raws. The simplest way to 
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eliminate this difference by choosing £ ee = 7r is not permitted by relations (|50|) . However, 
by the examination of Eqs.(f48|) and (|49| ) one can see that they are equivalent for £ eg = tt, 
£ee = igg = tt/2, £ se = and £, eg = (, gg = 7T / 2 , £ ge = £ ee = 0. Such choice does not contradict 
with Eqs.([45|), (|50|), it corresponds to £ 9 i = 0, £ g2 = ir/2, £ e i = n/2 , £ e2 = and £ 5 i = 
£32 = 7r/2, £ e i = ^ e 2 = 0. Inserting such values of £ Q/ 3 and 6>i = tt/4, 02 = — 7r/4 into Eqs.(^), 
(0) we find 

[ifx] = [M(0,0,7r/2,0)][i2(7r/4,-7r/4)]. (51) 

Otherwise, the matrix [i^] -1 is, by definition, given by Eqs.([H]). Inserting there Q\ = 7r/4, 
6*2 = 7r/2 and the values of £ a j given above we arrive to 

[H 2 ] = [R(—7r/ 4, — 7r/2)][M(0, — 7r/2, —tt/2, 0)]. (52) 

where we take into account that [i?(7r/4, 7r/2)] _1 = 7r/4, — 7r/2)] and 

[M(0, tt/2, tt/2,0)]- 1 = [M(0, -tt/2, -tt/2,0)]. Matrices [iy 1>2 ] are given explicitly by 
Eqs.(pn|). Obviously, that matrices given by Eqs.(^Tj), (^) are not the only ones which 
satisfy Eq.(^). However other possible local transformations will be similar to the ones 
given by matrices (|5l|) , (|52|) . 

ACKNOWLEDGEMENTS This work was supported by the European IST/FET 
project 'QUBITS' and by the European IHP network 'QUEST'. I. Protsenko is also grateful 
to Russian Foundation for Basic Research, grant 01-02-17330, for support. 



[1] A. Rauschenbeutel et al, Phys. Rev. Lett. 83, 5166 (1999). 

[2] G. Nogues et al, Nature 400, 239 (1999). 

[3] C. A. Sackett et al, Nature 404, 256 (2000). 

[4] A.Steane et al, Phys. Rev. A 62, 042305 (2000). 

[5] H. Walther, Proc. Roy. Soc. London A 454, 431 (1998). 

[6] N. A. Gershenfeld and I. L. Chuang, Science 275, 350 (1997). 

[7] A. Imamoglu et al, Phys. Rev. Lett. 83, 4204 (1999). 

[8] I. E. Protsenko, G. Reymond, N. Schlosser, and P.Grangier, Phys. Rev. A 65, 052301 (2002). 

[9] Ch.Roos et al, Phys. Rev. Lett. 83, 4713 (1999). 

[10] S. B. Zheng and G.C. Guo, Phys. Rev. Lett. 85, 2392 (2000); S. Osnaghi et al., Phys. Rev. 

Lett. 87, 037902 (2001). 

[11] C. J. Hood, and H. J. Kimble, Phys. Rev. A 64, 033804 (2001). 

[12] E. Knill, R. Laflamme and G.J. Milburn, Nature 409, 46 (2001). 



22 



[13] R. H. Dicke (unpublished); see also W. M. Itano et al, in Proceedings of the International 
Conference on Lasers '93 Lake Tahoe, 1993, edited by V. J. Corcoran and T. A. Goldman 
(STS Press, McLean, VA, 1994), pp.412 - 419. 

[14] C. Cabrillo, J. I. Cirac, P. Garcia-Fernandez, and P. Zoller, Phys. Rev. A 59, 1025 (1999). 

[15] N. Schlosser, G. Reymond, I. Protsenko, and P. Grangier, Nature 411, 1024 (2001). 

[16] J. H. Eberly, Am. J. Phys. 70, 276 (2002). 

[17] J. S. Bell, in Foundations of Quantum Mechanics, ed. by B. d'Espagnat (Academic, New York, 
1972). 

[18] A. Aspect, P. Grangier and G. Roger, Phys. Rev. Lett. 49, 91 (1982) 
[19] W. M. Itano et al, Phys. Rev. A 57, 4176 (1998). 
[20] G. Wilfred et al, Jap. Journ. Appl Phys, Part 1 40, 2100 (2001). 
[21] T. H. Oosterkamp et al, Nature, 395, 873 (1998). 



23 



